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The rogue wave solutions of the nonlinear Schrodinger equation (NLS) are tested in numerical 
simulations of weakly nonlinear and fully nonlinear hydrodynamic equations. Only the lowest order 
solutions from 1 to 5 are considered. A higher accuracy of wave propagation in space is reached using 
the modified NLS equation (MNLS) also known as the Dysthe equation. This numerical modelling 
allowed us to directly compare simulations with recent results of laboratory measurements in Q. In 
order to achieve even higher physical accuracy, we employed fully nonlinear simulations of potential 
Euler equations. These simulations provided us with basic characteristics of long time evolution of 
rational solutions of the NLS equation in the case of near breaking conditions. The analytic NLS 
solutions are found to describe the actual wave dynamics of steep waves reasonably well. 



INTRODUCTION 

Rogue waves in the ocean is an intriguing geophysical 
problem, which has attracted much interest in the recent 
years (see reviews 0^3] )• Although quasi-linear random 
wave superposition in principle can generate extremely 
high waves, it is the nonlinear effect of the Benjamin- 
Feir instability of deep-water surface waves which is now 
believed to be the main reason why the rogue wave oc- 
currence is so high in comparison with predictions of the 
basic Rayleigh theory Q . The intermediate stage of 
the instability is associated with generation of phase- 
correlated structures nonlinear wave groups, similar to 
that known within the framework of the integrable non- 
linear Schrodinger (NLS) equation 0-0- The NLS equa- 
tion is a very basic first approximation for the sea wave 
dynamics in the limit of small wave steepness and long 
wave modulations. It is amazing that the weakly nonlin- 
ear analytic solutions governed by the NLS turned out 
to be well reproducible in laboratory conditions even for 
the steep (almost breaking) waves [l|, Tol- T3|. 

One class of the NLS solutions has particular interest, 
since they are simple enough to be derived analytically, 
and at the same time they exhibit outstanding capability 
of wave enhancement. The Peregrine breather is the first 
wave solution among this hierarchy, which has become 
recognized as the general prototype of rogue waves on 
the background of ordinary waves [l3|, U | ■ It describes 



the growth of infinitesimal perturbation of a plane wave, 
which culminates in the amplification that is three times 
the initial wave amplitude; after a fleeting stage of ap- 
pearance, the rogue wave disappears and the plane wave 



is restored again (at least within the NLS equation frame- 
work). More general kinds of breather solutions were 
discovered in [15l - ll7| . 

The hierarchy of higher-order rational breather solu- 
tions has also been found 

(nHHl- The y describe a non- 
linear superposition of a group of rogue waves, which is 
localized both in time and space; they correspond to the 
case of degenerate eigenvalues of the associated scattering 
problem for the NLS. It was found (l9| that the maximum 
wave amplification described by the solutions is given by 
the formula 2N + I, where N is the order of the rational 
solution. The Peregrine breather corresponds to N = 1. 
Observation of rational solutions of the NLS equation in 
a laboratory has been reported in a series of publications 
[ToL UlL [22)]. Such waves have been named super rogue 
waves. Up to N = 5 solutions were successfully repro- 
duced in a wave flume Remarkably, the fifth-order 
solution exhibits the maximum wave amplification that 
is II times above the level of the background wave. 

The experiments in [l| demonstrated a good agreement 
between the analytic results given by the NLS solutions 
and the experimental observations. The agreement was 
achieved in both the dynamics of modulated waves and 
the observed shapes of the focused wave groups. The 
boundary condition for the wavemaker was specified in 
the form An (—d,t), where An is the analytic solution 
of the NLS equation, which corresponds to the A-order 
rational breather, and which attains the maximum am- 
plification at the position x = in the moment of time 
t = 0. Since the wave tank has limited length L, the 
following inequality \d\ < L must be satisfied to enable 
measuring of the wave profile at the location of the max- 
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imum amplitude (which corresponds to x — within the 
NLS framework). Indeed, the wave profiles measured at 
a distance « d downstream (near x = 0) were close to 
the analytic prototypes even when the focused wave was 
near the breaking point. This fact allows us to claim that 
the analytic rational breather solutions capture the non- 
linear effects well and can be considered as a good first 
approximation. 

At the same time, limitations of the experimental facil- 
ity (limited length of the flume, gauge accuracy, capillary 
and dissipation effects) resulted in the use of rather short 
distances of the wave evolution (d is small). The wave 
growth occurs in the scale of e~ 2 where e = ka is the wave 
steepness of the background wave. This is relatively long 
process in comparison with the propagation time along 
the tank. Since the distance d is small, the perturbation 
generated by the wavemaker is not infinitesimal but has 
finite amplitude. The wave enhancement from one side 
of the tank to the focusing point is only a part of the 
total amplification 2N + 1 predicted by the theoretical 
solution on a long distance. 

In the present study we analyse longer distances of 
wave propagation than those achievable in a water tank. 
This is done using numerical simulations rather than 
actual experiments. Namely, we solve partial differen- 
tial equations that are more accurate in describing wa- 
ter waves than NLS. Current results complement those 
presented in [lj. As in the experiment, only rational so- 
lutions of orders N = 1,2,3,4,5 are considered. The 
techniques and analytic expressions were given earlier in 
19 21]. The actual expressions are cumbersome and can- 
not be reproduced here. Two main questions posed in our 
previous work and addressed here are: 

1. How longer distances/times influence the nonlinear 
wave focusing in realistic conditions of steep waves 
(larger d and, correspondingly, r)? Do analytic so- 
lutions of the NLS provide adequate models for de- 
scription of nonlinear wave focusing of intense near 
breaking waves? 

2. What is the profile of the amplified wave near max- 
imum. Does the intensity of the background waves 
changes, when the focused waves start to break? 
The wave breaking onsets for breathers of differ- 
ent orders were estimated in the experiments [ljj. 
Namely, higher-order breathers N = 5 with ex- 
tremely small steepness of the background wave 
(ka > 0.01) were found to be breaking due to the 
nonlinear focusing. 

We examine the observations made in laboratory ex- 
periments [l[ and compare them firstly with the results 
of improved model for nonlinear wave modulations pre- 
sented in Section 2 and secondly with fully nonlinear sim- 
ulations given in Section 4. These techniques allowed us 
to make a better modelling of breathers appearing on top 



of surface gravity waves. We paid special attention, in 
Section 3, to the wave enhancements reached by the non- 
linear wave focusing in the weakly and strongly nonlinear 
approaches. Main outcomes of the study are summarised 
in the Conclusions. 



WAVE EVOLUTION WITHIN THE 
FRAMEWORK OF THE DYSTHE EQUATION 

The first-order approximation for deep-water surface 
gravity waves is the NLS in form 



— +C r — 

dt 9T dx 




A = 0. 



(1) 

Parameters k and uj are carrier wavenumber and cyclic 
frequency, which are linked according to the deep-water 
dispersion relation luq = gk$. The constant g is the 
gravity acceleration, while the group velocity is given by 
C gr = 2$- ■ The first approximation for the water surface 

elevation r)(x,t) is then given by 

r) = Re [A(x, t) exp(iu;o£ — ik x)]. 

More convenient for direct comparison with laboratory 
experiments are evolution equations that describe wave 
propagation in space. Using linear relations involving 
the group velocity, the NLS ([lj can be transformed to 
the following form 
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^L4| 2 A = 0. (2) 



From point of view of the asymptotic approach, Eq.([2]) 
is equally valid as ([T]). The latter form is predominantly 
used in optics. 

In this Section, we use the so called Dysthe model to 
simulate wave dynamics more realistically. This is the 
modified nonlinear Schrddinger (NLS) equation (MNLS), 
which takes into account higher-order terms responsible 
for nonlinear dispersion, full dispersion of linear waves 
and the effect of induced long-scale current [3, |24j . Du- 
ality of the form of the classic NLS equation (JXJ) , ©, 
allows us to write down either temporal and spatial ver- 
sions of the MNLS. For convenience, we use here the 
version of the MNLS equation, which describes the wave 
evolution in space 
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where: 
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0, for z < rj; 
0, for z — > c 



Just like in the laboratory runs, we use higher-order 
rational analytic solutions of the NLS A?q(x = —d,t), 
N = 1,2,3,4,5, to specify the boundary condition for 
the MNLS simulations starting at x = —d. Note that 
two different forms of the NLS equation (1) and (2) al- 
lows one the same solution to be written twofold. In 
order to start simulations of the spatial version of the 
MNLS equation ([3]), we use the corresponding solutions 
of the spatial version of the NLS equation @. The pro- 
gramming code of the Dysthe model used here has been 
verified earlier in the laboratory studies of intense water 
wave packets [25| . 

Similar to the classic NLS equation, the Dysthe model 
(3) governs the evolution of the complex wave amplitude 
A(x,t). Therefore it is a straightforward task to com- 
pare them with solutions of the NLS equation. As we 
used pseudospectral method to simulate the solutions of 
Dysthe equation, the time domain must have periodic 
boundary conditions. Although the breathers are de- 
fined on the infinite domain, they are localised and we 
used time intervals much wider than the characteristic 
localisation interval. 

Time series should be used to start the simulations 
with evolution in space, which is governed by equations 
Q . These time series can be easily calculated using ana- 
lytic expressions and are not given here. To give an idea, 
their wave profiles are similar to the space series shown 
in Fig. 4 except that independent variable is t. For con- 
sistency with experimental results, the wave amplitudes 
and periods as well as the associated wave lengths are 
chosen in accordance with those used in laboratory ex- 
periments All essential experimental parameters are 
given in Table HI 



Order of 
rational 
solution 


Steepness of 

the 
unperturbed 
wave, koAo 


Amplitude of 

the 
unperturbed 
wave Aq (m) 


Carrier 
wavelength 
A = If M 


1 


0.117 


0.01 


0.54 


2 


0.05 


0.003 


0.38 


3 


0.04 


0.002 


0.31 


4 


0.03 


0.003 


0.67 


5 


0.02 


0.002 


0.67 



TABLE I. Parameters for simulations of the MNLS equation. 

Having more freedom in conducting numerical experi- 
ments we used significantly longer distances of propaga- 
tion. Instead of the limiting length of d — 9 m in real 



experiment [l[ we used of up to 100 m in simulations. 
As can be seen from Table 1, the distance d — 100 m 
corresponds to the number of 150 - 320 carrier wave 
lengths. The results of these simulations are shown in 
Fig.l. Clearly, the wave amplification is the most fasci- 
nating feature of the breather evolution. The evolution of 
the maximum of the complex wave amplitude, i.e., max 
as a function of x is shown in Fig.l. This is done for 
five different orders of the rational breather, from N = 1 
to N = 5. Each simulation is started at the distance d to 
the left from the x = point with d — 9, 20, 30, 50 and 
100 m. The coding for each curve is shown in the left 
hand side insets. For comparison, the maximum ampli- 
tude of the exact breather solution is shown by the thin 
black curve. 

A close examination of Fig.l shows that despite a 
somewhat lower amplification observed in numerical sim- 
ulations, the curve shapes are quite similar to those given 
by the NLS solution (thin solid line in each panel). The 
difference becomes noticeable when N > 2 and at d val- 
ues higher than 50 m. This difference between the an- 
alytic NLS solution and the numerical solution of the 
Dysthe equation becomes qualitative in the case of the 
largest distance of d = 100 m. A substantial deviation in 
Fig.l can be observed for the orders 3, 4 and 5. 

Parasitic wave modulations grow along with the wave 
evolution. They are mainly located at the ends of the 
background wave train which is unavoidable feature of 
the technique and irrelevant to the breather itself. Ex- 
amples of these parasitic modulations can be seen in Fig. 5 
which is done for fully nonlinear simulations. The effect 
is very similar in the case of weakly nonlinear simulations 
of the Dysthe model. 

In order to eliminate this effect, the wave shifts close 
to the central peak should be carefully analysed when 
computing the maxima of \A\. Other unstable modes 
may also grow in the central area resulting in developing 
additional maxima in Figs.l. These maxima may appear 
at distances different from the expected location x — 
of the wave focusing. 

Another quite visible feature of Fig.l is that the case 
N = 1 is qualitatively different from N = 2, 5. There 
is a clear trend that at moderate values of d the amplifi- 
cation grows with d. Then the starting conditions which 
correspond to larger d result in a stronger wave amplifica- 
tion. When N = 1, even the longest distance d = 100 m 
results in strong wave amplification. This is not the case 
when N > 2. When N = 1, the maximum wave focusing 
actually occurs later than predicted by the NLS theory. 
Moreover, the delay increases with the distance d. On 
the other hand, there is no particular order in location 
of the focusing position in the cases of N > 2 although 
generally the maximum is closer to x = 0. 

Thus, the NLS equation provides good qualitative evo- 
lution of the point of maximum in the case of N = 1 but 
predicts its appearance earlier than it happens in the 
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FIG. 1. (Color online) Evolution of the maximum amplitude 
| -A | in numerical simulations of the Dysthe model. The insets 
show the distance d of propagation for each case before reach- 
ing the focusing point at x = 0. A thin black line in each panel 
corresponds to the analytic solution for the rational breathers 
of the NLS. The lower limit along the vertical axis shows the 
amplitude of the unperturbed background wave, Aq. The ex- 
perimental parameters that correspond to each case are given 
in Table 1. 



Dysthe model. The amplification values of w 3 for the 
Peregrine breather predicted in each model are in very 
good agreement. We should not forget that the starting 
boundary condition with small deviation from the con- 
stant amplitude background wave in simulations is the 



exact solution of the NLS equation. It may happen that 
corrections of this function may improve the agreement 
between the two models. 



WAVE AMPLIFICATION 

One of the parameters calculated in our study is the 
wave amplification \ defined as the ratio of the maximal 
amplitude A max to the amplitude of the background wave 
Aq, i.e. x = A max /Ao. In order to analyse this parameter 
in more detail, we plot it in Fig. 2. Namely, Fig. 2 displays 
the value of amplification reached in the simulations of 
the Dysthe model. This figure also contains the results 
of the strongly nonlinear simulations described below as 
well as experimental data. 

Generally, the curves in Fig.l have maxima at locations 
different from x = 0. Correspondingly, there are two sets 
of the same symbols in the plots - the larger and the 
smaller ones. The larger symbols show the maximum 
amplification reached within the simulation domain —d< 
x < 20 m, while the smaller symbols represent only the 
amplification reached in the close vicinity of the point 
x = 0. 

The circles in these plots show the ratio p\ of absolute 
maximum to the unperturbed envelope amplitude, i.e. 

pi = max \A\ I max \A-^(x — > oo, t)\ . 

The upward (red) and downward (blue) directed triangles 
show the wave enhancement in terms of surface elevation 
and depression, i.e., 

P2 = max(ry)/ max(|Ajv(a; — >• oo,£)|) 



and 



P3 



\rj\)/ min(|Ajv(:z: — > oo,i)| 



respectively. The surface elevation rj(x,t) is computed 
taking into account bound waves up to the third order. 
This approach is standard for the Dysthe theory, the de- 
tails can be found in 23[ 26 1. Horizontal lines in Fig. 2 
show the theoretical maxima 2N + 1 predicted by the 
NLS solutions. Generally, the complex wave amplitudes 
\A\ (circles) do not achieve the theoretical limit 2N + 1. 
However, quite remarkably, the contribution of bound 
wave components results in a greater amplification of 
wave crests max(^). These values can even exceed the 
theoretical limit 2N + 1 (upward directed triangles). 

When N = 1, the focusing is delayed. Then the max- 
imum which is measured at x = may be much smaller 
than after this point. Consequently, the smaller symbols 
show lesser amplification than the larger ones in the case 
N = 1. On the other hand, when N > 2, the maximum 
wave amplitudes in the entire domain — d < x < —20 m 
are similar to the ones near the point x = 0. 
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FIG. 2. (Color online) The relative values of wave maxima 
vs the distance d. Notations are as follows. Simulations of 
the Dysthe model: Wave amplification in terms of pi (circles 
o), p2 (upward directed triangles A) and pi (downward di- 
rected triangles The larger size symbols among circles 
and triangles are the maxima attained within the simulated 
distance (— d < x < 20 m). The smaller symbols of the same 
type denote the maxima in close proximity to the point x — 0, 
typically within one wavelength. Green asterisks (*) show the 
results of fully nonlinear simulations for near breaking condi- 
tions given in Section 4 and in the TableHTlwith the steepness 
values shown in bold. Crosses (x) correspond to the results 
of the experiment [l[ . The horizontal line in each panel at the 
level of 2N + 1 is the amplification provided by the analytic 
NLS solution. 
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FIG. 3. (Color online) Same as Fig.2, but modified amplifi- 
cations xi (circles o), (upward directed triangles A) and 
X3 (downward directed triangles v) are shown along the ver- 
tical axis. The upper limit of the vertical axis in the cases 
iV = 2, 3, 4, 5 corresponds to the amplification provided by 
the corresponding analytic NLS solution, 2N + 1. Green as- 
terisks ( ) show the results of fully nonlinear simulations for 
near breaking conditions. 



As can be seen from Fig.2, for N = 1 larger values of d 
result in larger amplifications. Vice versa, when N > 2, 
larger amplifications are observed at smaller d. The re- 
sults obtained in laboratory measurements [H are shown 
by crosses in Fig.2. In the experiment, d = 9 m and 
this point can be compared with small upward directed 
triangles at the same distance d. The numerical and ex- 



perimental points coincide in the case N = 1 while the 
experimental points are slightly higher when N = 2, 3, 4. 
When N = 5, the experimental point is lower than the 
one obtained in numerical simulations. This can be at- 
tributed to significantly smaller steepness, koAo = 0.01 
in the laboratory experiment for N = 5. We also note 
that the vertical scale in each panel does not start from 
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zero and the differences are smaller than it may seem at 
first glance. 

The results presented in Fig. 2 clearly demonstrate that 
the maximum surface elevation at x = obtained for the 
Dysthe model with d = — 9 m is very close to the value 
which is (27V + 1) times the amplitude of the background 
wave. These are given by the small upward directed tri- 
angles at d = —9 m. Close data were indeed observed 
in the experiments []]]. At the same time, the Dysthe 
model predicts the maximum wave amplification in the 
case of N = 1 to occur at a larger distance than the NLS 
theory provides. This feature still has to be confirmed in 
laboratory tests. 

The lower limit of the vertical axis in Fig.l is the back- 
ground wave amplitude Aq = An (x — > oo). When TV is 
large, i.e. AT = 3, 4, 5, the wave train is already amplified 
with respect to the background wave amplitude Aq at 
every x before the focusing point. The amplification is 
noticeable even at the largest distance d = 100 m consid- 
ered in simulations, see Fig.l. If we use the real meaning 
of 'amplification', i.e. its value relative to the starting 
wave amplitude, the actual numbers would be smaller. 

To take this correction into account, we calculated the 
amplifications with respect to the wave intensity at the 
starting point. Namely, we calculated the three following 
relative quantities 



Xi 



and 



X2 



X3 



max |^4| 


max \An (x = 


-d,t)\' 


max (77) 




max (j]n (x = 


-d,t)y 


min (77) 





min (r]N (x = —d, t)) 



instead of the those used in Fig. 2. The new characteristic 
parameters are presented in Fig. 3. Notations here are 
similar to those used in Fig. 2, i.e. the circles, upward and 
downward directed triangles have the same meaning. 

Clearly, if the initial amplitudes would coincide with 
the background one, Figures 2 and 3 would be identical. 
The data in Fig. 3 differ significantly from those in Fig. 2 
for short distances d and most visibly for large values 
of N. For the shortest focusing distance d = 9 m the 
modified wave amplification Xi is around 2. The maxi- 
mum wave amplification observed in all cases shown in 
Fig. 3 is a little bit larger than three. Even though the 
NLS solution with N > 2 describes the dynamics better 
for shorter distances d, the overall wave enhancement is 
higher when d is large. It is remarkable, that the first 
order solution exhibits systematically better agreement 
with the NLS exact solution when the distance d in- 
creases. The higher-order cases, N > 2 do not have a 
clear trend. 



ROGUE WAVES SIMULATED IN FULLY 
NONLINEAR EULER EQUATIONS 

Simulations according to the Dysthe model provide the 
most convenient data for comparison with the labora- 
tory experiments. However, the wave parameters given 
in Table 1 are very close to the breaking onset. The 
laboratory experiments [l| show that the steeper waves 
are indeed breaking. Therefore, using the approximate 
Dysthe model still can be inadequate. 

In this Section, we give the results of the numerical 
simulations of the Euler equations. We assume potential 
fluid motions. Most of the numerical runs are performed 
using the Higher-Order Spectral Method (HOSM). The 
details of the numerical approach can be found in [27j | . 
These simulations resolve up to 7-wave interactions, i.e. 
we use M = 6. Thus, practically, our numerical experi- 
ments are fully nonlinear. Selected runs of near-breaking 
focused waves were duplicated with the solver of the Eu- 
ler equations in conformal variables 28, 29[. The latter 
does not make any assumptions on wave steepness, and 
thus is a fully nonlinear approach. Comparing the two 
we can validate the results of the HOSM simulations. 
Consequently, we do not distinguish the results obtained 
using the two numerical approaches. 

The numerical codes have been developed previously 
in [30[. We are mainly interested in the strongly non- 
linear case when the focused wave is close to breaking. 
Near-breaking waves have the greatest enhancement rate 
31 1. The unstable mode which is the closest analog of 
the Peregrine breather with N = 1 has the amplification 
which actually exceeds the value of four. Remarkably, 
this is higher than the original amplification of 3. 

The steepness (&;o^4o) of the near-breaking NLS solu- 
tion was estimated on the basis of laboratory measure- 
ments in [![. This data are shown in the last column of 
Table 2. Based on the results of the previous Sections, 
one may expect that this threshold should be universal 
for the case N = 1. Indeed, the amplification does not 
depend very much on d. Moreover, it converges to a sin- 
gle value when d increases. 

For small distances d and N > 2, there is no universal 
steepness limit for wave breaking. This happens because 
the amplification factor and the wave steepness depend 
strongly on d as can be seen from Figs. 2 and 3. For large 
d, the breaking criterion should be close to the instability 
limit of the N = 1 case. 

The Euler equations govern the evolution in time. The 
approach is qualitatively different from solving the NLS 
@ or Dysthe ([3]) equations. The initial wave profile along 
x must be given. Posing the initial value problem re- 
sults in the fully nonlinear solution. We used several sets 
of initial conditions. We followed the laboratory exper- 
iments [H taking A]y(x,t = —L/C gr ) where L = 9 m 
is the effective length of the tank and C gr is the carrier 
group velocity. The focusing distance in terms of the 
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1 


16.76A w 33.52T 


0.090 


0.095 


0.12 


2 


23.87A w 47.74T 


0.060 


0.065 


0.06 


3 


28.65A » 53.30To 


0.040 


0.045 


0.05 


4 


14.32A w 28.64T 


0.030 


0.035 


0.04 


5 


28.65A « 53.30T,, 


0.025 


0.030 


0.02 



TABLE II. Parameters of the fully nonlinear simulations. 



wavelength, Ao and the wave period, To, is given in the 
second column of Table [TTJ The wavelengths A are given 
in Table 1. 

As we have to deal with periodic boundary conditions, 
the total numerical grid has been chosen much larger 
than the space interval we are interested in. This can be 
seen in Fig. 4. Various wave steepnesses fco^o have been 
considered for each wavelength. The intention was to 
find the steepness value when the breaking occurs. The 
latter manifests itself in abrupt instability of numerical 
simulations. 

The values of parameters when the wave evolution is ei- 
ther smooth or breaks down are given in the third column 
of Table HI1 The cases of survived waves with smooth evo- 
lution are highlighted in bold. The cases with breaking 
are shown in normal fonts. Thus, the breaking threshold 
parameters are in between these values. The breaking 
thresholds found in laboratory experiments [l[ are given 
in the last column of Table HB Comparison of the last two 
columns in Table lU shows that the breaking limit in the 
case N = 1 is much smaller in numerical simulations than 
in the corresponding laboratory measurements. One of 
the possible explanations is that in the experiments waves 
were not monitored at distances longer that d i.e. at 
x > 0. 

As can be seen from Figs. 1-3 of the previous Section, 
the strongest wave amplification in the case N = 1 oc- 
curs at x > 0. This observation agrees with the results 
of (3l|). Our fully numerical simulations were stopped 
either when the wave breaking had occurred or when the 
wave maximum had been reached. Slightly lower esti- 
mates of the steepness threshold for wavebreaking may 
occur for values N = 3,4. However, in general, the fully 
nonlinear simulations capture the wave breaking onset 
rather well. 

The second set of initial conditions corresponds to 
longer focusing times, when the initial conditions are cho- 
sen in the form An (x,t = — r), where r = 50To, lOOTo 
and 150Tp. The carrier wavenumbcr in all cases is set to 



fco = 20 rad/s, and the steepness is given in bold fonts in 
Table IE 

The amplifications observed in the fully nonlinear nu- 
merical simulations are shown in Figs. 2 and 3 with as- 
terisks. Finding the amplification value needs an appro- 
priate recalculation of the distance according to the dif- 
ference in the wavelengths of the carrier. In all cases the 
amplification values exceed the results of weakly nonlin- 
ear simulations of the Dysthe model but clearly follow 
the trends predicted by the Dysthe theory. They also ex- 
ceed the values obtained for near breaking waves in the 
laboratory experiment. 

The corrected maximum amplification of waves shown 
in Fig. 3 is around 3-4 times independent of the order of 
the rational solution N. Higher wave amplifications ob- 
served in fully nonlinear simulations in comparison to 
the weakly nonlinear simulations can be explained by 
the strongly nonlinear nature of the focused waves. The 
wavetrains at the time of the highest amplification are 
shown in Figs. 5 and 6. These are also presented for all 
five values of N. The initial conditions for these simula- 
tions are shown in Fig. 4. They are chosen at 100 wave 
periods prior the focusing point. The initial wave steep- 
ness in all cases is close to the breaking point. These are 
shown in bold in Table [TTJ 

The initial conditions shown by solid curves in Fig. 4 
correspond to the surface elevation in scaled coordinates. 
The dashed curves show the analytic NLS solutions for 
the wave envelope. The Stokes waves have vertical asym- 
metry making the wave crests slightly higher than the 
upper envelope, and the wave troughs shallower. This ef- 
fect is less pronounced for waves with smaller amplitudes. 
For larger N, the initial condition have more complicated 
shape. 

The bright curves in Fig. 5 show the surface eleva- 
tions at the moment when the maximum wave height is 
reached. The boundary of the hatched area clearly em- 
phasises the envelope of the wave. The hatching itself is 
done with several surface elevation curves slightly shifted 
in phase. The central region of the curves is zoomed and 
shown separately in Fig. 6. The main surface elevation 
curve here is contrasted in green while hatching is done 
by the black phase shifted lines just like in Fig. 5. Similar 
extreme wave trains were observed in simulations j3~0 ]. 
and very much similar solitary wave groups were mea- 
sured recently in laboratory conditions 
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To summarise, there is a set of wave packets with high 
amplitude central points that can be generated due to 
the nonlinear focusing. Just like in the NLS case, we 
can generate wave envelope structures with increasing 
number of nodes when increasing the values of N. Five 
lowest order wave profiles are shown in Fig. 5. These can 
be predicted based on NLS results although there is no 
complete correspondence with the NLS solutions. As ex- 
pected, the asymmetry is the main difference from the 
NLS case. Such asymmetry has been clearly observed 
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FIG. 4. (Color online) Initial conditions for the fully non- 
linear simulations. Surface elevation is shown by the solid 
curves. The envelopes described by the corresponding ana- 
lytic solutions of the NLS equation are shown by the dashed 
(red) curves. Each wavetrain is calculated 100 wave periods 
in advance of the wave focusing, r = 100 To. The background 
wave steepness is highlighted in bold in corresponding lines 
of Table M 



in the experiments In each case, the central region 
of the packet contains only a few wave oscillations. The 
central high amplitude part of the wave packet is an- 
other qualitative feature that in principle confirms the 
NLS prediction. Thus, the plots in Figs. 5 and 6 gen- 
erated in our numerical simulations qualitatively agree 
with the higher-order NLS rogue wave solutions as well 



FIG. 5. (Color online) Wavetrains at the point of the maxi- 
mum wave amplification obtained in the fully nonlinear sim- 
ulations with initial conditions shown in Fig.4. The bright 
curves show the surface elevation while the grey area fills the 
envelope. 



as with their observations in the laboratory 

An extensive and detailed study of modulated wave 
groups with high amplitude central peak is performed 
in [31| . The main conclusion made in 31 [ is that the 
nonlincarly focused waves, which are close to the wave 
breaking limit, have certain common features. The set of 
numerical simulations presented here also confirms this 
general conclusion. Firstly, the waves with maximum 
surface elevation shown by the thick green lines in Fig. 6 
have narrower crests than the Stokes wave. Secondly, the 
crest tips have the shape of the steepest Stokes wave (see 
details in [3l|). 



CONCLUSIONS 

The goal of our present work is to study the hierarchy 
of rogue wave solutions given by the exact solutions of the 
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FIG. 6. (Color online) The central region of the curves shown 
in Fig. 5 in a magnified scale. 



integrablc NLS in real life. We concentrated on solutions 
of this set with lowest order from TV = 1 to TV = 5. 
These solutions were recently observed experimentally in 
a water wave tank [![. The observations were found to 
agree reasonably well with the analytic solutions despite 
the limitations of the laboratory facility. This laboratory 
study motivated the research which is reported in the 
present paper. 

Namely, here we try to avoid the experimental limita- 
tions related to the short length of the wave tank and 
present simulations that describe the rogue waves with 
higher accuracy and in longer evolution than the exper- 
imental tank allows. Using two different approaches, we 



reproduced exact solutions of the integrable NLS. In one 
set of simulations we used Dysthc equations while the 
second technique is based on direct simulations of the 
Eulcr equations. 

Like in the experiment the starting conditions in 
our simulations are specified according to the exact NLS 
solutions. For higher accuracy, we took into account a 
few bound wave components. Further evolution then is 
compared with the NLS predictions. One of the main 
questions addressed by the simulations is how the growth 
of the rogue wave until the maximum amplitude is de- 
scribed in the two approaches. The main advantage of 
the strongly nonlinear simulations is the ability to de- 
scribe the evolution near the wave breaking limit. 

In general, we confirm that the wave evolution in new 
numerical simulations follows the exact breather solu- 
tions of the NLS equation reasonably well. The following 
features have been observed. 

1. There is a qualitative difference between the low- 
est TV = 1 and higher-order TV > 1 solutions. As 
expected, the dynamics described by the Peregrine 
breather (TV = 1) is much more robust than the 
dynamics prescribed by the solutions for TV > 1. 

2. The distance required for focusing in the simula- 
tions is longer than predicted by the analytic solu- 
tion in case TV = 1. For higher-order rogue waves 
TV > 1 the focusing distances in simulations and in 
theory are in good agreement. 

3. For the lowest order rogue wave TV = 1 the envelope 
amplification of three is well confirmed by our nu- 
merical simulations. On the other hand, solutions 
with TV > 1 provide noticeably smaller amplifica- 
tion than prescribed by the NLS theory. Even the 
wave evolution may differ if the focusing distance 
is too long. 

4. The effect of bound wave components is found to be 
very important for the process of focusing. Due to 
the contribution of bound waves, the wave crests 
are higher than the depth of wave troughs. Due 
to this difference, the enhancement of wave ampli- 
tudes is effectively larger than without the bound 
waves. In experiments, the amplification is usu- 
ally determined based on measurements of the wave 
crest amplitudes. 

5. Longer focusing distances lead to distortions of 
nonlinear wave focusing. As a result, the actual 
wave maxima for TV > 1 does not noticeably ex- 
ceed the value, which is observed in the case TV = 1 . 
The fully nonlinear simulations of nearly breaking 
waves show the maximum wave crest amplification 
up to about 4 times of the background. This esti- 
mate is in agreement with the results of numerical 
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simulations in [32[ and 3l[. Consequently, the an- 
alytic NLS solutions with TV > 1 are not optimal 
with respect to the maximal wave enhancement for 
use as initial conditions for nearly breaking waves. 
For waves with smaller amplitudes, the agreement 
with the NLS theory may be better. On the other 
hand, then the growth time or distance increases 
as « e -2 , where e is the wave steepness. 
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